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SI Dual-sugar mathematical model 

An ODE model of the GAL pathway was constructed based on extensive previous molecular characteri- 
zation of this system. This model was able to provide explanations for qualitative trends observed in the 
experimental data and insights into dynamics of the system. The model relied on the following simplifying 
assumptions: 

• No distinction was made between Gallp and Gal3p since these regulators have the same mechanism 
for activating the GAL genes by sequestration of Gal80p [1]. We assumed that these two species could 
be represented by a single species, Gallp. 

• Intracellular transport of glucose through the hexose transporters and galactose through Gal2p was 
not modeled. 

• A glucose repressor, R (for example Migl) transcriptionally represses GAL1 and GAL4 [2]. 

• Dimerization of Gal4p and Gal80p was not modeled [3, 4]. 

• We did not differentiate between nuclear and cytoplasmic partitioning of the GAL proteins [5, 6, 7]. 

The Hill coefficients for Gl (nci, nm), G4 (ura) and G80 (n 8 o) were estimated as 3, 2, 1 and 2 [8]. 
Based on these assumptions, the model that captures a set of molecular interactions for the wild-type (WT) 
GAL network in the presence of glucose and galactose is 
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where Gl, G4 and G80 represent Gallp, Gal4p and Gal80p. We assume that the the concentrations of the 
complexes, Gl-gal (Gl*), R-glu (R*), Gallp-Gal80p (C81) and Gal4p-Gal80p (C84) are at quasi-steady-state 
with respect to the dynamics of Gl, R, G4 and G80. Using this assumption, we can simplify the system of 
equations above to the following four ODEs 
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where Po = , Pi = -y c ^+ kr84 , P2 = ^+£ R , P3 - fc/ fl gal, Pi = Jois + k rg , p 5 = k m - k r8lPo , 

Pe = kf$ip 3 + k rS ip 0 p 3 , p 7 = pik rSi - fc/84 and p 8 = pik rSi - fc /84 . 

At steady-state, — = ^g-^ = = 0 and we could solve for the equilibrium concentrations of 

Gl, R and G4 in terms of G80. The equilibrium value of G80, G80 e was determined by computing the roots 
of a ninth order polynomial 

a 0 + a!G80 e H h a g G80 9 e = 0, 

where the coefficients, a,, are functions of the model parameters. The GAL80 feedback deletion model 
was constructed by replacing the Hill function of G4 by a constant acsas- The equilibrium points of this 
model were determined by computing the roots of a sixth order polynomial of G80 e . The stability of the 
equilibrium points was determined by computing the eigenvalues of the Jacobian matrix of the system of 
equations evaluated at each equilibrium point [9]. 

S2 Modification of 5 g and FoN-mid in a set of mutants 

Glucose repression is a complex process that modifies the activity of the GAL pathway at both the tran- 
scriptional and post-transcriptional levels [10]. Transcriptional repression of the GAL genes is mediated 
by a set of DNA binding proteins that recruit the global transcriptional repression complex Cyc8-Tupl, 
modifying the promoter state using multiple mechanisms including chromatin remodeling [11, 2, 12, 13, 14]. 
In response to galactose, the GAL promoters can escape a repressed state by the liberation of Gal4p from 
Gal80p and the inhibition of Cyc8-Tupl, which then initiates an ordered recruitment of the SAGA complex 
and TBP [15]. A PHD domain protein, Cti6p, has also been shown to play an important role in overcoming 
the transcriptional repression of the GAL1 promoter by Cyc8-Tupl [16]. 

In addition, the bifunctional glucose kinase, Hxk2p, also has been shown to play a critical role in glucose 
repression by regulating the nuclear to cytoplasmic ratio of Miglp by blocking its phosphorylation by the 
Snfl kinase [17, 18]. Hxk2p has a glucose-dependent nuclear localization and physically interacts with Miglp 
to regulate a set of target genes [19, 20, 21]. Furthermore, a previous study showed that a deletion of HXK2 
disrupts the hierarchical consumption order of glucose and galactose [22]. Interestingly, the recruitment of 
Cyc8-Tupl to the GAL1 promoter, as opposed to the changes in localization of Migl, has been shown to be 
the dominant regulatory interaction that forces this system into a repressed state [23] . 

To gain mechanistic insight into the role of glucose repression in generating the pre-emptive galactose 
pathway bimodal strategy (Fig. 1) , we characterized a set of mutants that modify the transcriptional state 
of the GAL promoters. For this experiment, we started with similar cell densities for each mutant (OD 60 o 
ss 0.3). Our results showed that a deletion of CTI6 significantly extended S g . Indeed, there were cells in 
the repressed subpopulation that never turned on the GAL genes over a period of approximately 14 hours 
(Fig. S4A-1). We also found that the fraction of cells in the ON state at the midpoint of the transient 
bimodal region (FoN-mid) was significantly lower in the CTI6A strain compared to WT (Fig. S4D). These 
data suggest that the timing of the repressed state's activation depends on CTI6 at the promoter- level. 

Interestingly, deletion of the dominant glucose kinase, Hxk2, significantly extended S g compared to WT 
(Fig. S4A-2). Similar to the CTI6A cells, the cells in the repressed subpopulation did not switch to the 
ON state during the course of the experiment. This mutant was particularly intriguing since it was the 
only modification to the system in which both FoN-mid and S g were augmented (Fig. S4C,D). Previous 
results have shown that cells deleted for Hxk2p consume glucose and galactose simultaneously [22]. Taken 
together, we hypothesize that Hxk2p is a crucial player that links the availability of glucose to the repression 
of galactose consumption and timing of GAL gene induction. 

We next characterized the dynamics of a strain in which Gal80p was regulated by an inducible TET pro- 
moter as opposed to its endogenous Gal4p-dependent promoter, therefore abolishing the feedback regulation 
of this protein (GAL80A fb) [8]. In the presence of 19.7 ng/ml aTc, the repressed subpopulation did not 
induce the GAL genes over the course of the experiment (Fig. S4A-3). These data demonstrate that a high 
and unregulated concentration of Gal80p relative to WT can block the repressed population from inducing 
the GAL genes at the appropriate time, highlighting the sensitivity of the network to perturbations in the 
regulation of this repressor. 
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We also identified three mutants that displayed smaller S g phenotypes compared to WT. The MIG1 
binding site is degenerate and other transcription factors can bind to the same sequence [24]. We hypoth- 
esized that a gene deletion of MIG1 could increase the probability that these alternative transcriptional 
regulators bind to the putative MIG1 sites and thus potentially alter the expression level of target genes. 
To directly reduce the promoter-level repression without significantly perturbing promoter occupancy, four 
point mutations were introduced into Migl (Migi4 m contains L490A, L493A, L496A and L498A) that have 
been shown to reduce the affinity of Migl to the Cyc8-Tupl complex [25, 23, 26]. This mutant exhibited a 
smaller 5 g and a larger FoN-mid compared to WT (Fig. S4B-1,C,D). 

We examined the role of the Cyc8-Tupl complex in the regulation of the transient bimodal dynamics. 
Since a deletion the Cyc8-Tupl complex has severe pleiotrophic effects, including flocculation, diploids were 
constructed that contain only one allele of CYC8 and TUP1 and therefore half the dosage of this complex. 
We compared the dynamics of this mutant to a wild type diploid and observed a significantly smaller S g and 
larger FoN-mid compared to wild type (Fig. S4B-2,C,D). 

To measure the contribution of Gal80p-Gal4p sequestration on the dynamics of the system, the affinity 
of Gal80p to Gal4p was reduced by mutating F856C or M861C of GAL4 [27]. The F856C mutation reduces 
the affinity of Gal80 to Gal4 more significantly than M861C. These two mutants exhibited a significantly 
smaller 6 g than WT (Fig. S4B-3,C). F856C had a smaller S g compared to M861C, suggesting that the 
binding affinity of the Gal80p-Gal4p complex is a crucial variable that controls the duration of bimodality. 
In addition, these two mutants exhibited a larger FoN-mid compared to WT (Fig. S4D). 

We also measured the growth rates of these strains to determine if the observed changes in the system's 
dynamics could be explained by alterations in the growth rates and thus the consumption rate of glucose. 
However, the changes in 5 g and FoN-mid were significantly greater than the variation in the growth rates 
(Fig. S4E). Therefore, these data strongly suggest that the observed modifications of the transient bimodality 
phenotype were dominated by the perturbations to the system's regulation. 

Taken together, these results suggest that the timing of the delayed activation response and the bias 
in the fraction of ON cells is regulated in a distributed way by many factors. Future investigation will 
reveal how these multiple players combine to generate the observed dynamic phenotypes. However, our data 
indicate that promoter-level regulation plays a critical role in controlling these phenotypes. 

53 Galactose does not benefit cells until glucose is consumed 

To assess the contribution of galactose to growth in the presence of glucose, we measured the extracellular 
concentration of glucose, galactose, growth and GAL gene expression dynamics of a strain lacking the 
endogenous GAL4 gene and expressing a DNA binding mutant of GAL4 (C14Y). This mutant renders 
Gal4p unable to activate transcription of the GAL genes [28] (Fig. S10A). Our data showed that this 
mutant consumed glucose at a rate equivalent to WT but was unable to consume galactose (Fig. S10B,C). 
Furthermore, the growth rates of the mutant and WT were very similar while the cells were consuming 
glucose (Fig. S10D). These data demonstrate that galactose does not benefit cells until glucose is fully 
consumed. 

54 Constitutive GAL gene expression reduces the glucose consump- 
tion rate 

Our data indicated that a high expression level of the GAL genes in a subpopulation of cells was not sufficient 
to initiate galactose metabolism (Fig. 2). To test whether altering GAL gene expression can modify the 
timing of galactose metabolism in the presence of glucose we used a strain in which the endogenous GALS 
gene was deleted and Gal3p was expressed from an inducible TET promoter regulated by an aTc-responsive 
transcription factor rtTA (GAL3A fb) [8, 29]. Using this strain, we induced the GAL system in advance of 
the wild type by administering 450 ng/ml aTc and monitored the dynamics of sugar consumption following 
induction of these two strains with 0.1% glucose and 0.1% galactose (Fig. S12A). 

The GAL3A fb strain exhibited an approximately 10% diminished rate of glucose consumption compared 
to wild type, indicating that constitutive GAL gene expression can have an inhibitory effect on glucose 
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metabolism (Fig. S12B). However, constitutive GAL gene expression was not sufficient to abolish the 
hierarchical order of sugar utilization. By contrast, galactose metabolism was accelerated in the mutant 
relative to WT (Fig. S12C). These combined effects reduced the metabolic delay, defined as the difference 
in the half-max of the galactose and glucose curves, by approximately 50% compared to WT (Fig. S12D). 

The GAL3A fb mutant displayed a reduced growth rate by up to 20% relative to WT (Fig. S12E,F). 
These data indicate that gratuitous expression of the GAL genes in environments where these genes are 
not required yields a fitness cost. This diminished growth and glucose consumption rate due to GAL gene 
overexpression provides insights into the benefits of tightly regulating the expression of the GAL genes in 
response to the availability of galactose in the environment. 

S5 Growth rates of the ON and OFF subpopulations 

For many dual-sugar inputs, the OFF and ON subpopulations were clearly distinguishable for several hours. 
We assumed that the transitions between the two states were insignificant since cellular switching between the 
OFF and ON states would yield intermediate fluorescence levels due to the gradual accumulation or dilution 
of fluorescent proteins. During this period, we computed the growth rates of the individual subpopulations 
as described in the Materials and Methods. Our results showed that the GAL OFF subpopulation grew on 
average 20% faster than the GAL ON subpopulation (Fig. S13). Since our data has shown that the glucose 
consumption rate can be reduced by GAL gene expression (Fig. SI 1 , SI 2), this effect could explain the 
diminished growth rate of the GAL OFF population in these conditions. 
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Strain name 

Wild-type (WT) P GA li Venus 

WT P GAL2 Venns 
WT P GAL7 Venus 
WT P GAL10 Venus 



Genotype 

MATaleu2::LEU2, trpl::TRPl-P GALJ Venus, ura3::URA3, ade2::ADE2, 

his3::HIS3 

MATaleu2::LEU2, trpl::TRPl-P Gj4L2 Venus, ura3::URA3, ade2::ADE2, 

his3::HIS3 

MATaleu2::LEU2, trpl::TRPl-P GAL7 Venus, ura3::URA3, ade2::ADE2, 
his3::HIS3 



MATa leu2::LEU2, 
ade2::ADE2, his3::HIS3 



trpl::TRPl-P GALJ0 Venus, ura3::URA3, 



WT P Gj4L2 Gal2- Venus, 

P G ALJomCherry 



MATa leu2, trpl, ura3::URA3-P Gj4 ijomCherry, 
his3::HIS3-P G/1 L2Gal2- Venus 



ade2::ADE2, 



Gal4-C14YP G/1L i 0 Venus 



MATa leu2::LEU2, trpl::TRPl-P GALJ0 Venus, ura3::URA3-P Gj4 ^Gal4- 
C14Y, ade2::ADE2, his3::HIS3, GAL4::HPH 



Gal4-ER P Gj4LJ0 Venus 



MATa leu2, trpl::TRPl-P G/1Li0 Venus, ura3 : : UR A3-P adhi G AL4D BD- 
ER-MSN2AD, ade2::ADE2, his3, GAL4::HPH 



HXK2A 



MATa leu2, trpl::TRPl-P G/1Li0 Venus, 
HXK2::NAT 



ura3, ade2::ADE2, his3, 



CTI6A 



MATa leu2, trpl::TRPl-P Gj4Li0 Venus, ura3, ade2::ADE2, his3, 
CTI6::NAT 



Migl 



4m 



MATa leu2::LEU2-P M / G iMIGl-L490A, L493A, L496A, L498A, 
trpl::TRPl-P Gj4Li0 Venus, ura3, ade2::ADE2, his3, MIG1::NAT 



Gal4 F856C 
Gal4 F861C 
WT diploid 

Cyc8-Tupl (+/-) diploid 
GAL3A fb* 
GAL80A fb* 
GAL2A fb* 



MATa leu2, trpl::TRPl-P GyiLJ0 Venus, ura3::URA3-P Gj4L ^GAL4- 
F856C, ade2::ADE2, his3, GAL4::HPH 

MATa leu2, trpl::TRPl-P Gj4LJ0 Venus, ura3::URA3-P Gj4L ^GAL4- 

F861C, ade2::ADE2, his3, GAL4::HPH 

MATa/a Ieu2/leu2, trpl::TRPl-P Gj4L2 Venus/trpl, ura3/ura3::URA3, 
ade2::ADE2/ade2, his3/his3 

MATa/a Ieu2/leu2, trpl::TRPl-P GAiJ Venus/trpl, ura3/ura3::URA3, 
ade2::ADE2/ade2, his3/his3, CYC8::NAT/CYC8, TUP1/TUP1::KAN 
MATa leu2::LEU2-P TBFm ^rtTA-M2, trpl::TRPl-P GALJ0 Venus, 
ura3::URA3-P TET GAL3, ade2::ADE2, his3, GAL3A::KAN 
MATa leu2::LEU2-P TBFm4 rtTA-M2, trpl::TRPl-P GALJ0 Venus, 

ura3 : : URA3-P t e tG AL80 , ade2::ADE2, his3, GAL80A::HPH 

MATa leu2::LEU2-P T£Fm4 rtTA-M2, trpl::TRPl-P GALJ0 Venus, 
ura3::URA3-P T BTGAL2, ade2::ADE2, his3, GAL2A::KAN 



Supplementary Table II: Strains used in this study. All strains were W303. * Strain described in [8]. 
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